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IN THE UNITED STATES ELECTED/DESIGNATED OFFICE 
OF THE UNITED STATES PATENT AND TRADEMARK OFFICE 
UNDER THE PATENT COOPERATION TREATY-CHAPTER H 

PRELIMINARY AMENDMENT 

APPLICANTS: Andreas Jurisch et al. DOCKET NO.: 112740-353 

SERIAL NO: GROUP ART UNIT: 

FILED: EXAMINER: 
INTERNATIONAL APPLICATION NO:: PCT/DE00/02434 
INTERNATIONAL FILING DATE 19 July 2000 

INVENTION: METHOD FOR DETERMINING THE AMPLITUDE AND 

PHASE ANGLE OF A MEASURING SIGNAL 
CORRESPONDING TO A CURRENT OR VOLTAGE OF 
AN ELECTRICAL POWER SUPPLY NETWORK 

Assistant Commissioner for Patents 
Washington, DC 20231 

PRELIMINARY AMENDMENT 

Sir: 

Please amend the above-identified International Application before entry into the 
National stage before the U.S. Patent and Trademark Office under 35 U.S.C. §371 as follows: 
In the Specification : 

Please replace the Specification of the present application, including the Abstract, 
with the following Substitute Specification: 

SPECIFICATION 
TITLE OF THE INVENTION 
METHOD FOR DETERMINING THE AMPLITUDE AND PHASE ANGLE OF A 
MEASURING SIGNAL CORRESPONDING TO A CURRENT OR VOLTAGE OF AN 
ELECTRICAL POWER SUPPLY NETWORK 
BACKGROUND OF THE INVENTION 
The present invention relates to a method for determining the amplitude and phase 
angle of a measuring signal corresponding to a current or a voltage on an electrical power 
supply network by using sampled values of the measuring signal, a model of the measuring 
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signal containing at least a sinusoidal component being used to calculate the model amplitude 
and phase parameters of the measuring signal with the sampled values by applying a recursive 
least-squares estimation method. 

A method of this type is described in an article by M.S. Sachdev and M. Nagpal "A 
recursive least squares error algorithm for power system relaying and measurement 
applications", IEEE Trans, on Power Delivery, Vol. 6, No. 3, July 1991. In the case of this 
known method, sampled values are formed from a measuring signal corresponding to a 
current or voltage of an electrical power supply network and, via a linear least-squares 
estimation method, by using a sinusoidal signal model that models the measuring signal, are 
used to form a complex vector which specifies the amplitude and phase angle of the 
measuring signal, hi this case, in a first step the real and imaginary part of the vector are in 
each case determined on their own. Then, in a second step, via a coordinate transformation, 
the polar coordinate representation of the complex vector, that is to say the magnitude and 
phase of the vector, can be determined from the real and imaginary part of the vector. 

In the case of the known method, it is assumed that the frequency of the measuring 
signal is known. If this is not so, or if the frequency changes, then a separate method is 
needed in order to determine the frequency of the measuring signal. Methods are known, for 
example, which measure the distance of the zero transitions of the measuring signal and, on 
the basis of measuring this period, determine the frequency of the measuring signal; see, for 
example, E. Schrufer (Editor): "Lexikon MeB- und Automatisierungstechnik" [Measurement 
and Automation Encyclopedia], VDI-Verlag, 1992, p. 204. Also known is a method for 
frequency measurement in which the measured signal to be examined is, in each case, filtered 
in parallel by a high-pass filter and an all-pass filter (German patent DE 42 1 1 946). The 
frequency of the measuring signal can be determined via the ratio between the amplitudes of 
the output signals from these two filters. 

The present invention is directed toward specifying a method with which all the 
significant variables of the measuring signal can be determined simultaneously and quickly. 
SUMMARY OF THE INVENTION 

In the case of a method of the type indicated at the beginning, this object is achieved, 
according to the present invention, in that use is made of a model of the measuring signal 
containing the sinusoidal component in accordance with the relationship y = A • sin(2 7ft + <p), 
y designating an instantaneous value of the model of the measuring signal, A the amplitude, f 



the frequency, 9 the phase angle and t the time. By using this model of the measuring signal 
and by using the sampled values, via a recursive nonlinear least-squares estimation method, 
the model frequency parameter of the measuring signal is also determined by the estimation 
together with the model amplitude parameter and the model phase angle parameter. 
5 Although it is known, from the book by H.J. Hermann, "Digitale Schutztechnik" 

[Digital Protection Engineering], 1997, p.p. 110-1 1 1, to use a recursive, nonlinear least- 
squares estimation method in protection engineering, the book does not contain any reference 
to the fact that by using such an estimation method, sampled values of a measuring signal can 
be used to determine the amplitude, phase angle and frequency of the measuring signal in a 
1 0 single measured-value processing process. 

A significant advantage of the method according to the present invention is, however, 
S precisely that the sampled values of the measuring signal are used to determine the frequency, 
R as well as the amplitude and phase angle, in a measured-value processing process. Therefore, 
!=& the amplitude, phase angle and frequency of the measuring signal are associated with the 
; 1;5 same point in time. 

!jy Furthermore, the patent DE 42 05 300 CI discloses a method with which the phase 

Q angle and the amplitude of a periodic signal can be determined via a phase-locked control 
jij loop (PLL (Phase-Locked Loop)). 

The use of the model for the measuring signal y = A • sin(2;rft + cp) leads to good 
Ilk) results if the measuring signal has a purely sinusoidal waveform. If there is a DC component 

present in the measuring signal, a model of the measuring signal in accordance with the 

relationship y = A • sin(2^ft + q>) + d is advantageously used, the summand d modeling the 

DC component of the measuring signal. 

If the measuring signal is a signal whose frequency changes over time, in the case of 
25 such measuring signals without a DC component, a model of the measuring signal in 

accordance with the relationship 

y = A-sm^2x]^(f {l} t')t + <pj 

advantageously can be used, and in the case of measuring signals with a DC component, a 
model of the measuring signal in accordance with the relationship 
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h'y + q) \ + d 



advantageously can be used, f^ designating the ith time derivative of the frequency and 
modeling a change in the frequency over time, and various orders of the time derivative of the 
frequency being taken into account by selecting the variable n. In these models, in an 
5 extension of the model specified above, y = A • sin(2^ft + cp), the frequency f is replaced by 

the expression ,=0 ; if only a Oth order derivative is taken into account, the integral 

expression becomes f^ or f. These expanded models make it possible to determine frequency 
changes £ l) over time, in addition to the variables amplitude A, phase angle cp and frequency f. 
, In a further advantageous embodiment of the method according to the present 

! 4=0 invention, the values of the amplitude A, the phase angle cp and the frequency f determined by 
yj the estimation method are output as results only when the estimation error is less than a 
vi smallest permitted estimation error. This has the advantage that, in particular, the values 

estimated at the start of the method and afflicted by large estimation errors are not output, and 
is therefore the large estimation errors cannot have any negative consequences for a user of the 
15 method. 

Additional features and advantages of the present invention are described in, and will 
Q be apparent from, the following Detailed Description of the Invention and the Figures. 

BRIEF DESCRIPTION OF THE FIGURES 
Figure 1 shows, by using a block diagram, the progress of an exemplary embodiment 
20 of the method according to the present invention. 

Figure 2 shows waveforms of a measuring signal, of an estimated signal y s determined 
for the measuring signal by using a model y = A • sin(2Mt + cp) and of the resulting error F. 

Figure 3 shows results of the method according to the present invention in the case of 
a purely sinusoidal measuring signal when use is made of an appropriate model of the 
25 measuring signal, separated in accordance with the time-based waveforms of amplitude, 
phase angle and frequency. 

Figure 4 shows results of the method according to the present invention in the case of 
a measuring signal subject to an offset, that is to say containing a DC component, separated in 
accordance with the time-based waveforms of amplitude, phase angle and frequency, when 
30 use is made of a model of the measuring signal that does not model any DC component. 
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Figure 5 shows results of the method according to the present invention in the case of 
a measuring signal subject to an offset, that is to say containing a DC component, separated in 
accordance with the time-based waveforms of amplitude, phase angle, frequency and DC 
component, when use is made of a model of the measuring signal that models a DC 
5 component. 

DETAILED DESCRIPTION OF THE INVENTION 
According to Figure 1, after being sampled in a sample-and-hold circuit 1 and after 
analog-digital conversion in an analog-digital converter 2, sampled values y m of a measuring 
signal u m are present on an input 3 of a data processing system 4. The data processing system 
10 4 contains a unit 5 for carrying out a recursive nonlinear least-squares estimation method, 
u At the beginning of the recursive estimation method, starting values SW for the 

%_ variables amplitude A, frequency f and phase angle q> are input into the unit 5 and are present 
:^ on the output 6 of the unit 5 as estimation output values © {o) . The estimation output values 

Uu 0 (Q) are led to a block 7, which contains a model of the measuring signal. 

15 In the block 7, the estimation output values @^ are used, taking into account the 

i=i | model of the measuring signal, to determine a starting signal value y s0 , which is output at the 
^ output of the block 7 and led to an input 8 of the unit 5. 

I U The unit 5 uses a (first) sampled value y m of the measuring signal u m and the starting 

signal value y s o to determine an estimation error F est im in accordance with the relationship (3) 

20 recited below. This estimation error will lie above a permitted smallest estimation error. 
Therefore, in the unit 5 and starting from the estimation output values @ (o) , in accordance 

with the rules for recursive nonlinear least-square estimations, new estimated values 0 (i) are 

formed in accordance with the relationship (1) below, are output at the output 6 of the unit 5 
and led to block 7. 

25 hi block 7, the new estimated values are used to determine an estimated signal 

value y s i, which is output at the output of block 7 and is led to the input 8 of the unit 5. The 
unit 5 uses the (first two) sampled values y m and the starting signal value y s0 and the 
estimated signal value y s i again to form an estimation error F est im in accordance with the 
relationship (3) recited below. This estimation error will generally also lie above the smallest 
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permitted estimation error. Therefore, in the unit 5 and starting from the estimated values 
on the preceding estimation, new estimated values @ (2) are again determined in 

accordance with the relationship (1). These steps are repeated until the estimation error F est im 
lies below a smallest permitted estimation error. 

In order to determine whether the estimation error F es tim lies below a smallest 
permitted estimation error, use is made of a block 9. Via an input 10, block 9 obtains from 
the unit 5 the sampled values y m of the measuring signal u m , the starting signal value y s0 and 
the estimated signal values y s i to y sk , which are buffered in the unit 5. When the estimation 
error F est i m lies below a smallest permitted estimation error, then a switching signal, which 
switches a switch 12 on, is output at an output 1 1 of the block 9. The switch 12 then 
forwards the estimated values present at the output 6 of the unit 5 to an output 13 of the 

e estin 

lZ from the estimation method and constitute the values for amplitude, frequency and phase 
ijj angle of the measuring signal, estimated with sufficient accuracy. After that, the estimation 
method continues to run as described above and permits continuous determination of the 



relationship (1). 

In relationship (1), @ w represents a vector which contains the estimated values of the 
variables to be determined, that is to say, here, the amplitude A, frequency f and phase angle 
cp, after k estimation steps; ® (k _ t) designates a vector which results in the case of the 

estimation after k-1 estimation steps. The matrix P (k) is a symmetrical precision matrix, as it 
is known, whose determination is presented further below using the relationship (2). The 

function ^<*> ^ contains the model of the measuring signal. The vector jq^ 



{i} 



contains the partial derivatives of the model {k) of the measuring signal with 
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respect to the parameters of the model of the measuring signal- that is to say, the derivatives 
with respect to amplitude, phase angle and frequency. The variable y (k > is the kth sampled 
value of the measuring signal. 

In order to determine the precision matrix P (k) , use is made of a relationship (2). 

i f ■ Lu t D ] (2: 

The factor X determines the exponential weighting of preceding sampled values of the 
measuring signal. At the beginning of the method, a starting value is used for the precision 
matrix P(k> 

The estimation error F estim after the kth estimation is determined via the relationship 

Here, N is the number to be taken into account of the sampled values y m of the 
measuring signal u m and of the estimated signal values y s , y mi is the ith sampled value of the 
measuring signal u m and y s i is the ith estimated signal value. 

The estimation error F es tim is determined by N sampled values y mi of the measuring 
signal u m and N estimated signal values y si being evaluated. An appropriate magnitude for N 
is the quotient of the (sampling frequency of the sample- and-hold circuit 1) / (estimated 
frequency), but N also can be selected to be larger. During the first estimation pass, the value 
used for N is the quotient of (sampling frequency of the sample-and-hold circuit 1) / (starting 
value SW of the frequency). When, at the beginning of the estimation method, N sampled 
values y mi of the measuring signal u m , or N estimated signal values y si are not yet present, then 
only the present values are used to determine the estimation error F es tim- 

Further details relating to carrying out recursive nonlinear least-squares estimation 
methods are described in the documents J. Wede; D. Werner: "EchtzeitprozeBmodelle auf der 
Basis von Parameterschatzverfahren" [Real-time process models based on parameter 
estimation methods], automation engineering series, Vol. 214, VEB Verlag Technik, Berlin, 
1985, p.p. 30-34, p.p. 44-50, 56, 57, and Hans- Joachim Herrmann: "Digitale Schutztechnik: 
Grundlagen, Software, Ausfuhrungsbeispiele" [Digital Protection Engineering: Principles, 
Software, Design Examples], VDE Verlag 1997, p.p. 104-113, and A. Jurisch: "Digitale 



ImpedanzmeBverfahren auf der Basis von Identifikationsmethoden" [Digital Impedance 
Measuring Methods Based on Identification Methods], Dissertation, Technische Hochschule 
Zittau, 1990, Annex 4.6.3. 

Figure 2 shows the time-based waveform of a sinusoidal measuring signal u m , whose 
5 amplitude, frequency and phase angle are to be determined by using the estimation method 
according to Figure 1 . Also shown over time, as the curve y s , are the estimated signal values 
y.i ... y S k which are output at the output of block 7. Also shown is the waveform of the values 
of the error F which results from the difference between the measuring signal u m and the 
estimated signal values y s i ... y S k, over time t. The waveforms shown result from an 
10 estimation method which is carried out by using the model y=A - sin(2^ + cp) for the 
U measuring signal. It can be seen that the waveform of the curve y s is the same as the 
y waveform of the measuring signal u m after about 20 ms; the error F then assumes very small 
W values. The estimation method, therefore, permits the correct determination of the amplitude, 
\j the frequency and the phase angle of the measuring angle after about one period of the 
~15 measuring signal u m . 

^ In Figure 3, taking into account the waveforms shown in Figure 2, the estimated 

l a values © ... linked with the estimated signal values y sJ ... y sk are plotted against the 

;jf time t, separated in accordance with the variables A, co and 27if and cp that determine them. In 
! u an upper graph, Figure 3 shows the waveform of the amplitude A determined by the 
20 estimation method according to Figure 1 , in a central graph the waveform of the variable 
co = 27tf that is proportional to the frequency f, and in a lower graph the waveform of the 
phase angle cp, in each case plotted against time. According to Figure 2, the result that the 
time-based waveforms have reached a steady state at their final value after about 20 ms with 
the smallest permitted estimation error, and indicate the result of the estimation method. 
25 Figure 4 shows, in a manner of representation which is the same as Figure 3, the time- 

based waveforms of amplitude A, frequency f and phase angle (p determined in a further 
estimation by using the model of the measuring signal y = A • sm(2nft + cp) if use is made of 
sampled values of a measuring signal which contains a DC component. The waveforms of 
frequency and phase angle also contain large time-based fluctuations for times greater than 
30 20 ms. These fluctuations point to an erroneous estimation and are a reference to the fact that 
the selected model is not beneficial for a measuring signal containing a DC component. 



Figure 5 shows the time-based waveforms of amplitude A, frequency f, phase angle <p 
and DC component d determined in an estimation using the model y = A • sin(2/rft + cp) + d 
for the measuring signal, if use is made of sampled values of a measuring signal which 
contains a DC component. By using the model just mentioned of the measuring signal, 
5 fluctuations in the time-based waveforms of frequency and phase angle are avoided; 

consequently, an estimation which is better with regard to the estimation error is carried out. 

Although the present invention has been described with reference to specific 
embodiment, those of skill in the art will recognize that changes may be made thereto without 
departing from the spirit and scope of the invention as set forth in the hereafter appended 
10 claims. 

ABSTRACT OF THE DISCLOSURE 
Q A method for determining the amplitude and phase angle of a measuring signal 

yj corresponding to a current or a voltage on an electrical power supply network by using 
C sampled values of the measuring signal, a model of the measuring signal containing at least a 

! =f 5 sinusoidal component being used to calculate the amplitude and phase of the measuring signal 

W 

with the sampled values by applying a recursive least-squares estimation method. 
|T In order to be able to determine the frequency of the measuring signal together with 

|]J the amplitude and the phase angle, use is made of a model of the measuring signal in 
Q accordance with the relationship y = A • sin(2;rft + cp), and by using this model and the 
20 sampled values (y m ), the determination of the frequency of the measuring signal (u m ) as well 

is carried out via a recursive nonlinear least-squares estimation method. By expanding the 

signal model, measuring signals with a DC component and with frequencies that change over 

time also can be investigated. 

In the claims: 

25 On page 12, cancel line 1 and substitute the following left-hand justified heading 

therefore: 

CLAIMS 

Please cancel claims 1-5, without prejudice, and substitute the following claims 
therefore. 

30 6. A method for determining an amplitude and phase angle of a measuring signal 

corresponding to a current or a voltage on an electrical power supply network by using 
sampled values of the measuring signal, the method comprising the steps of: 

9 



using a model of the measuring signal containing a sinusoidal component in 
accordance with the relationship y = A • sin(2^ + (p), y designating an instantaneous value of 
the model of the measuring signal, A designating the amplitude, f designating the frequency, 
cp designating the phase angle and t the time; and 
5 using both the model of the measuring signal and the sampled values, via a 

recursive nonlinear least-squares estimation method, to determine the model frequency 
parameter of the measuring signal by the estimation together with the model amplitude 
parameter and the model phase angle parameter. 

10 7. A method for determining an amplitude and phase angle of a measuring signal 

;_=, as claimed in claim 6, wherein use is made of a model for the measuring signal in accordance 
y with the relationship y = A ■ sm(2xft + cp) + d, d modeling a DC component of the measuring 

Q 

W signal. 

!=* 

%j 

'15 8. A method for determining an amplitude and phase angle of a measuring signal 

]l _ as claimed in claim 6, wherein use is made of a model for the measuring signal in accordance 



y = A-sm\2nY J (f {,) t')t + (p\ 

V /=0 J f(i) A 



m with the relationship ^ '' =0 ' , f*° designating the ith order time 

rj derivative of the frequency and modeling a change in the frequency over time, and various 

orders of the time derivative of the frequency being taken into account by selecting the 
20 variable n. 

9. A method for determining an amplitude and phase angle of a measuring signal 
as claimed in claim 7, wherein use is made of a model of the measuring signal in accordance 

y = A ■ sinl 2n J] (f m t ' > + cp + d 
with the relationship ^ ,=0 ' , f° } designating the ith time derivative of 

25 the frequency and modeling a change in the frequency over time, and various orders of the 

time derivative of the frequency being taken into account by selecting the variable n. 

10. A method for determining an amplitude and phase angle of a measuring signal 
as claimed in claim 6, wherein the values of the amplitude, the phase angle and the frequency 
determined by the estimation method are output as resulting values only when an estimation 

30 error is less than a smallest permitted estimation error. 
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REMARKS 

The present amendment makes editorial changes and corrects typographical errors in 
the specification, which includes the Abstract, in order to conform the specification to the 
requirements of United States Patent Practice. No new matter is added thereby. 



claims by the current amendment. The attached page is captioned " Versions w ith Markings 
to Show Changes Made." 

hi addition, the present amendment cancels original claims 1-5 in favor of new claims 
6-10. Claims 6-10 have been presented solely because the revisions by crossing out 
10 underlining which would have been necessary in claims 1-5 in order to present those claims 

in accordance with preferred United States Patent Practice would have been too extensive, 
M and thus would have been too burdensome. The present amendment is intended for 
Pj clarification purposes only and not for substantial reasons related to patentability pursuant to 
uy 35 U.S.C. §§101, 102, 103 or 1 12. Indeed, the cancellation of claims 1-5 does not constitute 
*§5 an intent on the part of the Applicants to surrender any of the subject matter of claims 1-5. 
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Attached hereto is a marked-up version of the changes made to the specification and 



Early consideration on the merits is respectfully requested. 



William E. Vaughan [/ 
Bell, Boyd & Lloyd LLC 
P.O.Box 1135 

Chicago, Illinois 60690-1135 
(312) 807-4292 
Attorneys for Applicants 




(Reg. No. 39,056) 
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VERSIONS WITH MARKINGS TO SHOW CHANGES MADE 



SPECIFICATION 
TITLE OF THE INVENTION 
5 METHOD FOR DETERMINING THE AMPLITUDE AND PHASE ANGLE OF A 

MEASURING SIGNAL CORRESPONDING TO A CURRENT OR VOLTAGE OF AN 
ELECTRICAL POWER SUPPLY NETWORK 
BACKGROUND OF THE INVENTION 

Description 

10 Method for Determining the Amplitude and Phase Angel of a Measuring Signal 
Corresponding to a Current or Voltage of an Electrical Power Supply Network 

The present invention relates to a method for determining the amplitude and phase 
angle of a measuring signal corresponding to a current or a voltage on an electrical power 
supply network by using sampled values of the measuring signal, a model of the measuring 

1 5 signal containing at least a sinusoidal component being used to calculate the model amplitude 
and phase parameters of the measuring signal with the sampled values by applying a recursive 
least-squares estimation method. 

A method of this type is described in an article by M.S. Sachdev and M. Nagpal "A 
recursive least squares error algorithm for power system relaying and measurement 

20 applications", IEEE Trans, on Power Delivery, Vol. 6, No. 3, July 1991 . In the case of this 
known method, sampled values are formed from a measuring signal corresponding to a 
current or voltage of an electrical power supply network and, by means of via a linear least- 
squares estimation method, by using a sinusoidal signal model that models the measuring 
signal, are used to form a complex vector^ which specifies the amplitude and phase angle of 

25 the measuring signal. In this case, in a first step the real and imaginary part of the vector are,, 
in each case,, determined on their own. Then, in a second step, by means of via a coordinate 
transformation, the polar coordinate representation of the complex vector, that is to say the 
magnitude and phase of the vector, can be determined from the real and imaginary part of the 
vector. 

30 In the case of the known method, it is assumed that the frequency of the measuring 

signal is known. If this is not so, or if the frequency changes, then a separate method is 
needed in order to determine the frequency of the measuring signal. Methods are known 



the*e, for example, which measure the distance of the zero transitions of the measuring signal 
and, on the basis of measuring this period, determine the frequency of the measuring signals- 
see, for example, E. Schrufer (Editor): "Lexikon MeB- und Automatisierungstechnik" 
[Measurement and Automation Encyclopedia], VDI-Verlag, 1992, p. 204. Also known is a 
5 method for frequency measurement in which the measured signal to be examined is s in each 
case., filtered in parallel by a high-pass filter and an all-pass filter (German patent 
DE 42 1 1 946). The frequency of the measuring signal can be determined via the ratio 
between the amplitudes of the output signals from these two filters. 

The present invention is har.ed on tho object of directed toward specifying a method 

10 with which all the significant variables of the measuring signal can be determined 
simultaneously and quickly. 

M SUMMARY OF THE INVENTION 

Q In the case of a method of the type indicated at the beginning, this object is achieved, 

I": according to the present invention, in that use is made of a model of the measuring signal 
15 containing the sinusoidal component in accordance with the relationship y = A ■ sin(2^ + <p), 
yj y designating an instantaneous value of the model of the measuring signal, A the amplitude, f 
□ the frequency, <p the phase angle and t the time^b v. By using this model of the measuring 
:r! signal and by using the sampled values, by means of via a recursive nonlinear least-squares 
It) estimation method, the model frequency parameter of the measuring signal is also determined 
ilo by the estimation together with the model amplitude parameter and the model phase angle 
parameter. 

Although it is known, from the book by H.J. Hermann, "Digitale Schutztechnik" 
[Digital Protection Engineering], 1997, p.p. 110-1 1 1, to use a recursive, nonlinear least- 
squares estimation method in protection engineering, the book does not contain any reference 
25 to the fact that by using such an estimation method, sampled values of a measuring signal can 
be used to determine the amplitude, phase angle and frequency of the measuring signal in a 
single measured-value processing process. 

A significant advantage of the method according to the present invention is, however, 
precisely that the sampled values of the measuring signal are used to determine the frequency, 
30 as well as the amplitude and phase angle, in a measured- value processing process^-and 

therefore . Therefore, the amplitude, phase angle and frequency of the measuring signal are 
associated with the same point in time. 



Furthermore, the patent DE 42 05 300 CI discloses a method with which the phase 
angle and the amplitude of a periodic signal can be determined via a phase-locked control 
loop (PLL (Phase-Locked Loop)). 

The use of the model for the measuring signal y = A • sin(2^ft + <p) leads to good 
results if the measuring signal has a purely sinusoidal waveform. If there is a DC component 
present in the measuring signal, a model of the measuring signal in accordance with the 
relationship y = A • sin(2^ft + cp) + d is advantageously used, the summand d modeling the 
DC component of the measuring signal. 

If the measuring signal is a signal whose frequency changes over time, in the case of 
such measuring signals without a DC component, a model of the measuring signal in 
accordance with the relationship 



can advantageously can be used, and in the case of measuring signals with a DC component, a 
model of the measuring signal in accordance with the relationship 

y-- 



©an advantageously can be used, f® designating the ith time derivative of the frequency and 
modeling a change in the frequency over time, and various orders of the time derivative of the 
frequency being taken into account by selecting the variable n. In these models, in an 
extension of the model specified above, y = A • sin(2/rft + cp), the frequency f is replaced by 

the expression 1=0 ; if only a 0th order derivative is taken into account, the integral 

expression becomes f (0) or f. These expanded models make it possible to determine frequency 
changes i® over time, in addition to the variables amplitude A, phase angle cp and frequency f. 

In a further advantageous refinem e nt embodiment of the method according to the 
present invention, the values of the amplitude A, the phase angle cp and the frequency f 
determined by the estimation method are output as results only when the estimation error is 
less than a smallest permitted estimation error. This has the advantage that, in particular, the 
values estimated at the start of the method and afflicted by large estimation errors are not 

14 



output, and therefore the large estimation errors cannot have any negative consequences for a 
user of the method. 

Additional features and advantages of the present invention are described in, and will 
be apparent from, the following Detailed Description of the Invention an d the Figures. 
BRIEF DESCRIPTION OF THE FIGURES 

In order to explain the invention further, Eigure 1 shows, by using a block diagram, 
the progress of an exemplary embodiment of the method according to the present invention^ 

Figure 2 shows waveforms of a measuring signal, of an estimated signal y s determined 
for the measuring signal by using a model y = A • sin(2^ft + q>) and of the resulting error F ?i 

Figure 3 shows results of the method according to the present invention in the case of 
a purely sinusoidal measuring signal when use is made of an appropriate model of the 
measuring signal, separated in accordance with the time-based waveforms of amplitude, 
phase angle and frequency^. 

Figure 4 shows results of the method according to the present invention in the case of 
a measuring signal subject to an offset, that is to say containing a DC component, separated in 
accordance with the time-based waveforms of amplitude, phase angle and frequency, when 
use is made of a model of the measuring signal that does not model any DC componemyand. 

Figure 5 shows results of the method according to the present invention in the case of 
a measuring signal subject to an offset, that is to say containing a DC component, separated in 
accordance with the time-based waveforms of amplitude, phase angle, frequency and DC 
component, when use is made of a model of the measuring signal that models a DC 
component. 

DETAILED DESCRIPTION OF THE INVENTION 
According to Figure 1, after being sampled in a sample-and-hold circuit 1 and after 
analog-digital conversion in an analog-digital converter 2, sampled values y m of a measuring 
signal u m are present on an input 3 of a data processing system 4. The data processing system 
4 contains a unit 5 for carrying out a recursive nonlinear least-squares estimation method. 

At the beginning of the recursive estimation method, starting values SW for the 
variables amplitude A, frequency f and phase angle q> are input into the unit 5 and are present 
on the output 6 of the unit 5 as estimation output values . The estimation output values 

(S\ are led to a block 7, which contains a model of the measuring signal. 
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In the block 7, the estimation output values ©^ are used, taking into account the 
model of the measuring signal, to determine a starting signal value y s0 , which is output at the 
output of the block 7 and led to an input 8 of the unit 5. 

The unit 5 uses a (first) sampled value y m of the measuring signal u m and the starting 
signal value y s o to determine an estimation error F est i m in accordance with the relationship (3) 
recited below. This estimation error will lie above a permitted smallest estimation error. 
Therefore, in the unit 5 and starting from the estimation output values © (q) , in accordance 

with the rules for recursive nonlinear least-square estimations, new estimated values are 
formed in accordance with the relationship (1) below, are output at the output 6 of the unit 5 
and led to block 7. 

In block 7, the new estimated values © are used to determine an estimated signal 

value y s i, which is output at the output of block 7 and is led to the input 8 of the unit 5. The 
unit 5 uses the (first two) sampled values y m and the starting signal value y s0 and the 
estimated signal value y s i again to form an estimation error F est im in accordance with the 
relationship (3) recited below. This estimation error will generally also lie above the smallest 
permitted estimation error. Therefore, in the unit 5 and starting from the estimated values 
© i on the preceding estimation, new estimated values ©^ are again determined in 

accordance with the relationship (1). These steps are repeated until the estimation error F estim 
lies below a smallest permitted estimation error. 

In order to determine whether the estimation error F es tim lies below a smallest 
permitted estimation error, use is made of a block 9. Via an input 10, block 9 obtains from 
the unit 5 the sampled values y m of the measuring signal u m , the starting signal value y s o and 
the estimated signal values y s i to y S k, which are buffered in the unit 5. When the estimation 
error F es tim lies below a smallest permitted estimation error, then a switching signal, which 
switches a switch 12 on, is output at an output 1 1 of the block 9. The switch 12 then 
forwards the estimated values present at the output 6 of the unit 5 to an output 13 of the 

data processing system 4. The estimated values ©^ are then output as the resulting values 

from the estimation method and constitute the values for amplitude, frequency and phase 
angle of the measuring signal, estimated with sufficient accuracy. After that, the estimation 
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method continues to run as described above and permits continuous determination of the 
estimated values . 

The unit 5 determines the estimated values @ w in accordance with the following 
relationship (1). 

® w =©(*-,) +£(*)£(*, ta) -% ( *,^(*-oJl 

hi relationship (1), @ represents a vector which contains the estimated values of the 
variables to be determined, that is to say, here, the amplitude A, frequency f and phase angle 
q>, after k estimation steps; © (fc _ i) designates a vector which results in the case of the 
estimation after k-1 estimation steps. The matrix P (k) is a symmetrical precision matrix, as it 
is known, whose determination is presented further below using the relationship (2). The 

function ^ (i > ' — ^ contains the model of the measuring signal. The vector j m 

contains the partial derivatives of the model ^ of the measuring signal with 

respect to the parameters of the model of the measuring signal; that is to say, the derivatives 
with respect to amplitude, phase angle and frequency. The variable y w is the kth sampled 
value of the measuring signal. 

In order to determine the precision matrix P^), use is made of a relationship (2). 

i r • loo t P 1 

P t . , = — p.. - = Y £(k-n 

The factor X determines the exponential weighting of preceding sampled values of the 
measuring signal. At the beginning of the method, a starting value is used for the precision 
matrix P (k ). 

The estimation error F est im after the kth estimation is determined by means of via the 
relationship 
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Here, N is the number to be taken into account of the sampled values y m of the 
measuring signal u m and of the estimated signal values y s , y mi is the ith sampled value of the 
measuring signal u m and y si is the ith estimated signal value. 
5 The estimation error F es ti m is determined by N sampled values y mi of the measuring 

signal u m and N estimated signal values y si being evaluated. An appropriate magnitude for N 
is the quotient of the (sampling frequency of the sample-and-hold circuit 1) / (estimated 
frequency), but N ean also can be selected to be larger. During the first estimation pass, the 
value used for N is the quotient of (sampling frequency of the sample-and-hold circuit 1) / 
CjO (starting value SW of the frequency). When, at the beginning of the estimation method, N 
i]j sampled values y mi of the measuring signal u m , or N estimated signal values y S i are not yet 
H present, then only the present values are used to determine the estimation error F est im- 
UJ Further details relating to carrying out recursive nonlinear least-squares estimation 

r methods are described in the documents J. Wede; D. Werner: "EchtzeitprozeBmodelle auf der 
hi"5 Basis von Parameterschatzverfahren" [Real-time process models based on parameter 
I « estimation methods], automation engineering series, Vol. 214, VEB Verlag Technik, Berlin, 
Q 1985, p.p. 30-34, p.p. 44-50, 56, 57, and Hans-Joachim Herrmann: "Digitale Schutztechnik: 
! y Grundlagen, Software, Ausfuhrungsbeispiele" [Digital Protection Engineering: Principles, 

Software, Design Examples], VDE Verlag 1997, p.p. 104-1 13, and A. Jurisch: "Digitale 
20 ImpedanzmeBverfahren auf der Basis von Identifikationsmethoden" [Digital Impedance 

Measuring Methods Based on Identification Methods], Dissertation, Technische Hochschule 
Zittau, 1990, Annex 4.6.3. 

Figure 2 shows the time-based waveform of a sinusoidal measuring signal u m , whose 
amplitude, frequency and phase angle are to be determined by using the estimation method 
25 according to Figure 1. Also shown over time, as the curve y s , are the estimated signal values 
y s i — y S k which are output at the output of block 7. Also shown is the waveform of the values 
of the error F which results from the difference between the measuring signal u m and the 
estimated signal values y s i ... y S k, over time t. The waveforms shown result from an 
estimation method which is carried out by using the model y=A - sin(2^? + (p) for the 
30 measuring signal. It can be seen that the waveform of the curve y s is the same as the 



waveform of the measuring signal u m after about 20 ms; the error F then assumes very small 
values. The estimation method^ therefore,, permits the correct determination of the amplitude, 
the frequency and the phase angle of the measuring angle after about one period of the 
measuring signal u m . 

5 In Figure 3, taking into account the waveforms shown in Figure 2, the estimated 

values © (i) ... © (t) linked with the estimated signal values y s i ...y sk are plotted against the 

time t, separated in accordance with the variables A, co and 2%f and cp that determine them. In 
an upper graph, Figure 3 shows the waveform of the amplitude A determined by the 
estimation method according to Figure 1, in a central graph the waveform of the variable 
10 © = 27cf that is proportional to the frequency f, and in a lower graph the waveform of the 
u phase angle cp, in each case plotted against time. According to Figure 2, the result that the 

° time-based waveforms have reached a steady state at their final value after about 20 ms with 

O 

UJ the smallest permitted estimation error, and indicate the result of the estimation method, 
ly Figure 4 shows, in a manner of representation which is the same as Figure 3, the time- 

fi 5 based waveforms of amplitude A, frequency f and phase angle cp determined in a further 
L estimation by using the model of the measuring signal y=A- sm(2 nft + <p) if use is made of 
M sampled values of a measuring signal which contains a DC component. The waveforms of 
frequency and phase angle also contain large time-based fluctuations for times greater than 
9 20 ms. These fluctuations point to an erroneous estimation and are a reference to the fact that 
20 the selected model is not beneficial for a measuring signal containing a DC component. 

Figure 5 shows the time-based waveforms of amplitude A, frequency f, phase angle cp 
and DC component d determined in an estimation using the model y = A • sin(2^ft + cp) + d 
for the measuring signal, if use is made of sampled values of a measuring signal which 
contains a DC component. By using the model just mentioned of the measuring signal, 
25 fluctuations in the time-based waveforms of frequency and phase angle are avoided; 

consequently, an estimation which is better with regard to the estimation error is carried out. 

Although the present invention has been described with reference to specific 
embodiment, those of skill in the art will recognize that changes maybe made thereto without 
departing from the spirit and scope of the invention as set forth in the h ereafter appended 
30 claims. 
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Method for determining the amplitude and phase angle of 
a measuring signal corresponding to a current or 
5 voltage of an electrical power supply network 

The invention relates to a method for determining the 
amplitude and phase angle of a measuring signal 
corresponding to a current or a voltage on an 

10 electrical power supply network by using sampled values 
of the measuring signal, a model of the measuring 
signal containing at least a sinusoidal component being 
used to calculate the model amplitude and phase 
parameters of the measuring signal with the sampled 

15 values by applying a recursive least-squares estimation 
method. 



|« A method of this type is described in an article by 

M.S. Sachdev and M. Nagpal "A recursive least squares 
□ 2 0 error algorithm for power system relaying and 
Ef s measurement applications", IEEE Trans. on Power 

ft} Delivery, Vol. 6, No. 3, July 1991. In the case of this 

known method, sampled values are formed from a 
measuring signal corresponding to a current or voltage 

2 5 of an electrical power supply network and, by means of 

a linear least-squares estimation method, by using a 
sinusoidal signal model that models the measuring 
signal, are used to form a complex vector, which 
specifies the amplitude and phase angle of the 
30 measuring signal. In this case, in a first step the 
real and imaginary part of the vector are in each case 
determined on their own. Then, in a second step, by 
means of a coordinate transformation, the polar 
coordinate representation of the complex vector, that 

3 5 is to say the magnitude and phase of the vector, can be 

determined from the real and imaginary part of the 
vector . 



GR 99 P 4121 

- 2 - 

In the case of the known method, it is assumed that the 
frequency of the measuring signal is known. If this is 
not so, or if the frequency changes, then a separate 
method is needed in order to determine the frequency of 
the measuring signal. Methods are known there, for 
example, which measure the distance of the zero 
transitions of the measuring signal and, on the basis 
of measuring this period, determine the frequency of 
the measuring signal, see, for example, E. Schriifer 
(Editor) : "Lexikon MeS- und Automatisierungstechnik" 
[Measurement and Automation Encyclopedia] , VDI-Verlag, 
1992, p. 2 04. Also known is a method for frequency 
measurement in which the measured signal to be examined 
is in each case filtered in parallel by a high-pass 
filter and an all-pass filter (German patent 
DE 42 11 946) . The frequency of the measuring signal 
can be determined via the ratio between the amplitudes 
of the output signals from these two filters. 

The invention is based on the object of specifying a 
method with which all the significant variables of the 
measuring signal can be determined simultaneously and 
quickly. 

In the case of a method of the type indicated at the 
beginning, this object is achieved, according to the 
invention, in that use is made of a model of the 
measuring signal containing the sinusoidal component in 
accordance with the relationship y = A ■ sin(2?rf~t + q>) , 
y designating an instantaneous value of the model of 
the measuring signal, A the amplitude, f the frequency, 
9 the phase angle and t the time; by using this model 
of the measuring signal and by using the sampled 
values, by means of a recursive nonlinear least - squares 
estimation method, the model frequency parameter of the 
measuring signal is also determined by the estimation 
together with the model amplitude parameter and the 
model phase angle parameter. 
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Although it is known, from the book by H.-J. Hermann, 
"Digitale Schutztechnik" [Digital Protection 

Engineering], 1997, pp. 110-111, to use a recursive, 
nonlinear least- squares estimation method in protection 
5 engineering, the book does not contain any reference to 
the fact that by using such an estimation method, 
sampled values of a measuring signal can be used to 
determine the amplitude, phase angle and frequency of 
the measuring signal in a single measured-value 
10 processing process. 

A significant advantage of the method according to the 
invention is, however, precisely that the sampled 
values of the measuring signal are used to determine 
15 the frequency, as well as the amplitude and phase 
angle, in a measured-value processing process, and 
therefore the amplitude, phase angle and frequency of 
the measuring signal are associated with the same point 
in time. 

20 

Furthermore, the patent DE 42 05 3 00 CI discloses a 
method with which the phase angle and the amplitude of 
a periodic signal can be determined via a phase -locked 
control loop (PLL (Phase-Locked Loop) ) . 

25 

The use of the model for the measuring signal 
y = A • sin(2;rft + cp) leads to good results if the 
measuring signal has a purely sinusoidal waveform. If 
there is a DC component present in the measuring 
3 0 signal, a model of the measuring signal in accordance 
with the relationship y = A • sin(2^ft + cp) + d is 
advantageously used, the summand d modeling the DC 
component of the measuring signal. 
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If the measuring signal is a signal whose frequency 



signals without a DC component, a model of the 
measuring signal in accordance with the relationship 



can advantageously be used, and in the case of 
measuring signals with a DC component, a model of the 
measuring signal in accordance with the relationship 



10 can advantageously be used, f (i) designating the ith 
time derivative of the frequency and modeling a change 
in the frequency over time, and various orders of the 
time derivative of the frequency being taken into 
account by selecting the variable n. In these models, 

15 in an extension of the model specified above, 
y = A • sin(2;rft + q>) , the ^ frequency 

f is replaced by the expression ^(/ ( 'V) '■ if only 
a 0th order derivative is taken ,=0 into 
account, the integral expression becomes f (0) or f. 

20 These expanded models make it possible to determine 
frequency changes f (i) over time, in addition to the 
variables amplitude A, phase angle cp and frequency f . 

In a further advantageous refinement of the method 
25 according to the invention, the values of the amplitude 
A, the phase angle cp and the frequency f determined by 
the estimation method are output as results only when 
the estimation error is less than a smallest permitted 
estimation error. This has the advantage that, in 
30 particular, the values estimated at the start of the 
method and afflicted by large estimation errors are not 
output, and therefore the large 



changes over time, 



in the case of such measuring 



5 
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estimation errors cannot have any negative consequences 
for a user of the method. 

In order to explain the invention further, 
5 figure 1 shows, by using a block diagram, the progress 
of an exemplary embodiment of the method according to 
the invention, 

figure 2 shows waveforms of a measuring signal, of an 
estimated signal y s determined for the measuring signal 
10 by using a model y = A • sin(27zft + <p) and of the 
resulting error F, 

figure 3 shows results of the method according to the 
invention in the case of a purely sinusoidal measuring 
signal when use is made of an appropriate model of the 
15 measuring signal, separated in accordance with the 
time-based waveforms of amplitude, phase angle and 
frequency, 

figure 4 shows results of the method according to the 
invention in the case of a measuring signal subject to 

2 0 an offset, that is to say containing a DC component, 

separated in accordance with the time-based waveforms 
of amplitude, phase angle and frequency, when use is 
made of a model of the measuring signal that does not 
model any DC component , and 
25 figure 5 shows results of the method according to the 
invention in the case of a measuring signal subject to 
an offset, that is to say containing a DC component, 
separated in accordance with the time -based waveforms 
of amplitude, phase angle, frequency and DC component, 

3 0 when use is made of a model of the measuring signal 

that models a DC component. 

According to figure 1, after being sampled in a sample - 
and-hold circuit 1 and after analog-digital conversion 
35 in an analog-digital converter 2, sampled values y m of 
a measuring signal u m are present on an input 3 of a 
data processing system 4 . The data processing system 4 



GR 99 P 4121 



- 6 - 



contains a unit 5 for carrying out a recursive 
nonlinear least -squares estimation method. 

At the beginning of the recursive estimation method, 
5 starting values SW for the variables amplitude A, 
frequency f and phase angle cp are input into the unit 5 
and are present on the output 6 of the unit 5 as 
estimation output values 0 (q) . The estimation output 

values 0 (q) are led to a block 7, which contains a model 

10 of the measuring signal. 

In the block 7, the estimation output values @ (o) are 

used, taking into account the model of the measuring 
signal, to determine a starting signal value y s0 , which 
15 is output at the output of the block 7 and led to an 
input 8 of the unit 5. 

The unit 5 uses a (first) sampled value y m of the 
measuring signal u m and the starting signal value y s0 to 

2 0 determine an estimation error F est im in accordance with 

the relationship (3) recited below. This estimation 
error will lie above a permitted smallest estimation 
error. Therefore, in the unit 5 and starting from the 
estimation output values , in accordance with the 

25 rules for recursive nonlinear least-square estimations, 
new estimated values are formed in accordance with 

the relationship (1) below, are output at the output 6 
of the unit 5 and led to block 7. 

3 0 In block 7, the new estimated values © (j) are used to 

determine an estimated signal value y s i, which is output 
at the output of block 7 and is led to the input 8 of 
the unit 5. The unit 5 uses the (first two) sampled 
values y ra and the starting signal value y s0 and the 
35 estimated signal value y si again to form an estimation 
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error F es tim in accordance with the relationship (3) 
recited below. This estimation error will generally 
also lie above the smallest permitted estimation error. 
Therefore, in the unit 5 and starting from the 
estimated values on the preceding estimation, new 

estimated values ©_ (2) are again determined in accordance 
with the relationship (1) . These steps are repeated 
until the estimation error F est im lies below a smallest 
permitted estimation error. 

In order to determine whether the estimation error F es tim 
lies below a smallest permitted estimation error, use 
is made of a block 9. Via an input 10, block 9 obtains 
from the unit 5 the sampled values y m of the measuring 
signal u m , the starting signal value y s0 and the 
estimated signal values y sX to y sk , which are buffered 
in the unit 5. When the estimation error F estiin lies 
below a smallest permitted estimation error, then a 
switching signal, which switches a switch 12 on, is 
output at an output 11 of the block 9. The switch 12 
then forwards the estimated values present at the 

output 6 of the unit 5 to an output 13 of the data 
processing system 4. The estimated values are then 

output as the resulting values from the estimation 
method and constitute the values for amplitude, 
frequency and phase angle of the measuring signal, 
estimated with sufficient accuracy. After that, the 
estimation method continues to run as described above 
and permits continuous determination of the estimated 
values © (jt) . 

The unit 5 determines the estimated values @ w in 
accordance with the following relationship (1) . 
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In relationship (1) , ®^ represents a vector which 

contains the estimated values of the variables to be 
determined, that is to say, here, the amplitude A, 
frequency f and phase angle q>, after k estimation 
5 steps; ®(k i) designates a vector which results in the 

case of the estimation after k-1 estimation steps. The 
matrix P< k ) is a symmetrical precision matrix, as it is 
known, whose determination is presented further 
below using the relationship (2) . The function ^(2<*)'— 

10 contains the model of the measuring signal. The 

vector y_(k) contains the partial derivatives of the 
model ^l^(t) » — (*-!) j °^ t ^ Le measuring signal with 

respect to the parameters of the model of 

the measuring signal, that is to say the derivatives 

15 with respect to amplitude, phase angle and frequency. 
The variable y^ is the kth sampled value of the 
measuring signal. 

In order to determine the precision matrix P (Jc) , use is 
2 0 made of a relationship (2) . 

The factor X determines the exponential weighting of 

2 5 preceding sampled values of the measuring signal. At 

the beginning of the method, a starting value is used 
for the precision matrix P(k) . 

The estimation error F est im after the kth estimation is 

3 0 determined by means of the relationship 



(3) 



GR 99 P 4121 



- 9 - 



Here, N is the number to be taken into account of the 
sampled values y m of the measuring signal u m and of the 
estimated signal values y s , y m i is the ith sampled value 
of the measuring signal u m and y s i is the ith estimated 
5 signal value . 

The estimation error F est im is determined by N sampled 
values y m i of the measuring signal u m and N estimated 
signal values y si being evaluated. An appropriate 
10 magnitude for N is the quotient of the (sampling 
frequency of the sample-and-hold circuit 1) / 
(estimated frequency) , but N can also be selected to be 

□ larger. During the first estimation pass, the value 

used for N is the quotient of (sampling frequency of 

\2 15 the sample-and-hold circuit 1) / (starting value SW of 

the frequency) . When, at the beginning of the 

Li 

jVj estimation method, N sampled values y m ± of the measuring 

signal u m , or N estimated signal values y S i are not yet 
;~: present, then only the present values are used to 

ill 2 0 determine the estimation error F es tim- 

Rj 
Q 

iil Further details relating to carrying out recursive 

nonlinear least-squares estimation methods are 
described in the documents J. Wede ; D. Werner: 

25 "Echt zeitproze&modelle auf der Basis von 

Parameterschatzverf ahren" [Real-time process models 
based on parameter estimation methods] , automation 
engineering series, Vol. 214, VEB Verlag Technik, 
Berlin, 1985, p.p. 30-34, pp. 44-50, 56, 57, and Hans- 

3 0 Joachim Herrmann: "Digitale Schut ztechnik : Grundlagen, 
Software, Ausf uhrungsbei spiel e" [Digital Protection 
Engineering: Principles, Software, Design Examples], 
VDE Verlag 1997, pp. 104-113, and A. Jurisch: "Digitale 
ImpedanzmeSverf ahren auf der Basis von 

3 5 Identif ikationsmethoden" [Digital Impedance Measuring 
Methods Based on Identification Methods] , Dissertation, 
Technische Hochschule Zittau, 1990, Annex 4.6.3. 
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Figure 2 shows the time-based waveform of a sinusoidal 
measuring signal u m , whose amplitude, frequency and 
phase angle 
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are to be determined by using the estimation method 
according to figure 1. Also shown over time, as the 
curve Ysi are the estimated signal values y s i • • . Ysk 
which are output at the output of block 7. Also shown 
5 is the waveform of the values of the error F which 
results from the difference between the measuring 
signal u m and the estimated signal values y s i . . . y S k, 
over time t. The waveforms shown result from an 
estimation method which is carried out by using the 

10 model y = A • sin(2^ft + (p) for the measuring signal. It 
can be seen that the waveform of the curve y s is the 
same as the waveform of the measuring signal u m after 
about 2 0 ms; the error F then assumes very small 
values. The estimation method therefore permits the 

15 correct determination of the amplitude, the frequency 
and the phase angle of the measuring angle after about 
one period of the measuring signal u m . 

In figure 3, taking into account the waveforms shown in 
20 figure 2, the estimated values ... 0 (i) linked with 

the estimated signal values y si . . . y S k are plotted 
against the time t, separated in accordance with the 
variables A, co and 27tf and q> that determine them. In an 
upper graph, figure 3 shows the waveform of the 

25 amplitude A determined by the estimation method 
according to figure 1, in a central graph the waveform 
of the variable eo = 27if that is proportional to the 
frequency f, and in a lower graph the waveform of the 
phase angle cp, in each case plotted against time. 

30 According to figure 2, the result that the time-based 
waveforms have reached a steady state at their final 
value after about 2 0 ms with the smallest permitted 
estimation error, and indicate the result of the 
estimation method. 

35 

Figure 4 shows, in a manner of representation which is 
the same as figure 3, the time-based waveforms 
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of amplitude A, frequency f and phase angle <p 
determined in a further estimation by using the model 
of the measuring signal y = A • sin(2;rft + q>) if use is 
made of sampled values of a measuring signal which 
5 contains a DC component. The waveforms of frequency and 
phase angle also contain large time-based fluctuations 
for times greater than 20 ms . These fluctuations point 
to an erroneous estimation and are a reference to the 
fact that the selected model is not beneficial for a 
10 measuring signal containing a DC component. 

Figure 5 shows the time-based waveforms of amplitude A, 
frequency f, phase angle <p and DC component d 
determined in an estimation using the model 

15 y = A • sin(2;rft + cp) + d for the measuring signal, if 
use is made of sampled values of a measuring signal 
which contains a DC component. By using the model just 
mentioned of the measuring signal, fluctuations in the 
time -based waveforms of frequency and phase angle are 

20 avoided; consequently, an estimation which is better 
with regard to the estimation error is carried out. 



GR 99 P 4121 



Patent claims 

1. A method for determining the amplitude and phase 
angle of a measuring signal corresponding to a 
current or a voltage on an electrical power supply 
network by using sampled values of the measuring 
signal, a model of the measuring signal containing 
at least a sinusoidal component being used to 
calculate the model amplitude and phase parameters 
of the measuring signal with the sampled values by 
applying a recursive least -squares estimation 
method, characterized in that 

use is made of a model of the measuring signal 
containing the sinusoidal component in 
accordance with the relationship 

y = A ■ sin(2^ft + q>) , y designating an 
instantaneous value of the model of the 
measuring signal, A the amplitude, f the 
frequency, 9 the phase angle and t the time, 
and 

by using this model of the measuring signal and 
by using the sampled values, by means of a 
recursive nonlinear least -squares estimation 
method, the model frequency parameter (f) of 
the measuring signal is also determined by the 
estimation together with the model amplitude 
parameter (A) and the model phase angle 
parameter (q>) . 

2. The method as claimed in claim 1, characterized in 
that use is made of a model for the measuring 
signal in accordance with the relationship 
y = A ■ sin(2^ft + cp) + d, d modeling the DC 
component of the measuring signal . 



The method as claimed in claim 1, characterized in 
that use is made of a model for the measuring 
signal in accordance with the relationship 
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= ,4- sin 2tf£(y c, V> + < 



f (i) designating the ith 
order time derivative of the frequency and 
modeling a change in the frequency over time, and 
various orders of the time derivative of the 
frequency being taken into account by selecting 
the variable n. 

The method as claimed in claim 2, characterized in 

that use is made of a model of the measuring 

signal in accordance with the relationship 

y = A-sm\2nJ^(f {i) t i )t + (p) + d 

v '= c ' , f u> designating the ith time 

derivative of the frequency and modeling a change 

in the frequency over time, and various orders of 

the time derivative of the frequency being taken 

into account by selecting the variable n. 

The method as claimed in one of claims 1 to 4 , 
characterized in that the values of the amplitude 
(A) , the phase angle (cp) and the frequency (f ) 
determined by the estimation method are output as 
resulting values only when the estimation error is 
less than a smallest permitted estimation error. 



GR 99 P 4121 



Abstract 

Method for determining the amplitude and phase angle of 
a measuring signal corresponding to a current or 
voltage of an electrical power supply network 

The invention relates to a method for determining the 
amplitude and phase angle of a measuring signal 
corresponding to a current or a voltage on an 
electrical power supply network by using sampled values 
of the measuring signal, a model of the measuring 
signal containing at least a sinusoidal component being 
used to calculate the amplitude and phase of the 
measuring signal with the sampled values by applying a 
recursive least - squares estimation method. 
In order to be able to determine the frequency of the 
measuring signal together with the amplitude and the 
phase angle, use is made of a model of the measuring 
signal in accordance with the relationship 
y = A ■ sin(2^ft + (p) , and by using this model and the 
sampled values (y m ) , the determination of the frequency 
of the measuring signal (uj as well is carried out by 
means of a recursive nonlinear least-squares estimation 
method. By expanding the signal model, measuring 
signals with a DC component and with frequencies that 
change over time can also be investigated. 



FIG 1 



ABSTRACT OF THE DISCLOSURE 

Method for determining the amplitude and phase angle of a measuring signal corresponding 
to a current or voltage of an electrical power supply network 
5 The invention relates to a A method for determining the amplitude and phase angle of 

a measuring signal corresponding to a current or a voltage on an electrical power supply 
network by using sampled values of the measuring signal, a model of the measuring signal 
containing at least a sinusoidal component being used to calculate the amplitude and phase of 
the measuring signal with the sampled values by applying a recursive least-squares estimation 
10 method. In order to be able to determine the frequency of the measuring signal together with 
y. the amplitude and the phase angle, use is made of a model of the measuring signal in 
'S accordance with the relationship y = A ■ sin(2;rft + cp), and by using this model and the 
W sampled values (y m ), the determination of the frequency of the measuring signal (u m ) as well 
'"•-J is carried out by means of via a recursive nonlinear least-squares estimation method. By 
|vi5 expanding the signal model, measuring signals with a DC component and with frequencies 
^ that change over time ean also can be investigated. 
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Declaration and Power of Attorney For Patent Application 
Erklarung Fur Patentanmeldungen Mit Vollmacht 

German Language Declaration 



Als nachstehend benannter Erfinder erklare ich hiermit 
an Eides Statt: 



dass mein Wohnsitz, meine Postanschrift, und meine 
Staatsangehorigkeit den im Nachstehenden nach 
meinem Namen aufgefuhrten Angaben entsprechen, 



dass ich, nach bestem Wissen der ursprungliche, erste 
und alleinige Erfinder (fails nachstehend nur ein Name 
angegeben ist) oder ein ursprunglicher, erster und 
Miterfinder (falls nachstehend mehrere Namen 
aufgefuhrt sind) des Gegenstandes bin, fur den dieser 
Antrag gestellt wird und fur den ein Patent beantragt 
wird fur die Erfindung mit dem Titel: 

Verfahren zum Ermitteln von Amplitude 
und Phasenwinkel eines einemStrom 
oder einer Spannunq eines elektrischen 
Eneraieversorqunqsnetzes 
** "entsprechenden Messsiqnals 

deren Beschreibung 

(zutreffendes ankreuzen) 
f~l hier beigefugt ist. 

am 19.07.2000 als 
PCT internationale Anmeldung PCT/DE00/02434 

PCT Anmeldungsnummer 

eingereicht wurde und am 

abgeandert wurde (falls tatsachlich abgeandert). 



Ich bestatige hiermit, dass ich den Inhalt der obigen 
Patentanmeldung einschliesslich der Anspruche 
durchgesehen und verstanden habe, die eventuell 
durch einen Zusatzantrag wie oben erwahnt abgean- 
dert wurde. 



Ich erkenne meine Pflicht zur Offenbarung irgendwel- 
cher Informationen, die fur die Priifung der vorliegen- 
den Anmeldung in Einklang mit Absatz 37, Bundes- 
gesetzbuch, Paragraph 1.56(a) von Wichtigkeit sind, 
an. 



Ich beanspruche hiermit auslandische Prioritatsvorteile 
gemass Abschnitt 35 der Zivilprozessordnung der 
Vereinigten Staaten, Paragraph 119 aller unten ange- 
gebenen Auslandsanmeldungen fur ein Patent oder 
eine Erfindersurkunde, und habe auch alle Auslands- 
anmeldungen fur ein Patent oder eine Erfindersurkun- 
de nachstehend gekennzeichnet, die ein Anmelde- 
datum haben, das vor dem Anmeldedatum der 
Anmeldung liegt, fur die Prioritat beansprucht wird. 



As a below named inventor, I hereby declare that: 



My residence, post office address and citizenship are 
as stated below next to my name, 



I believe I am the original, first and sole inventor (if only 
one name is listed below) or an original, first and joint 
inventor (if plural names are listed below) of the 
subject matter which is claimed and for which a patent 
is sought on the invention entitled 



Method for determining the amplitude 
and angle of phase of a measuring 
signal corresponding to the current 
or voltage of an electric power 
distribution system 

the specification of which 

(check one) 

□ is attached hereto. 

was filed on 19.07.2000 as 
PCT international application PCT/DE00/02434 

PCT Application No. 

and was amended on 

(if applicable) 



I hereby state that I have reviewed and understand the 
contents of the above identified specification, including 
the claims as amended by any amendment referred to 
above. 



I acknowledge the duty to disclose information which is 
material to the examination of this application in 
accordance with Title 37, Code of Federal Regulations, 
§1 .56(a). 



I hereby claim foreign priority benefits under Title 35, 
United States Code, §119 of any foreign application(s) 
for patent or inventor's certificate listed below and have 
also identified below any foreign application for patent 
or inventor's certificate having a filing date before that 
of the application on which priority is claimed: 
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Patent and Trademark Office-U.S. DEPARTMENT OF COMMERCE 



Prior foreign appplications 
Prioritat beansprucht 



19934055.2 DE 
(Number) (Country) 
(Nummer) (Land) 



(Number) 
(Nummer) 



(Number) 
(Nummer) 



(Country) 
(Land) 



(Country) 
(Land) 



German Language Declaration 



19.07.1999 

(Day Month Year Filed) 

(Tag Monat Jahr eingereicht) 



(Day Month Year Filed) 
(Tag Monat Jahr eingereicht) 



(Day Month Year Filed) 
(Tag Monat Jahr eingereicht) 



□ 
No 
Nein 



□ 
Yes 



□ 
Yes 



Ich beanspruche hiermit gemass Absatz 35 der Zivil- 
prozessordnung der Vereinigten Staaten, Paragraph 
120, den Vorzug aller unten aufgefiihrten Anmel- 
dungen und falls der Gegenstand aus jedem Anspruch 
dieser Anmeldung nicht in einer fruheren 
amerikanischen Patentanmeldung laut dem ersten 
Paragraphen des Absatzes 35 der ZivilprozeRordnung 
der Vereinigten Staaten, Paragraph 122 offenbart ist, 
erkenne ich gemass Absatz 37, Bundesgesetzbuch, 
Paragraph 1 .56(a) meine Pflicht zur Offenbarung von 
Informationen an, die zwischen dem Anmeldedatum 
der fruheren Anmeldung und dem nationalen Oder PCT 
internationalen Anmeldedatum dieser Anmeldung 
bekannt geworden sind. 



I hereby claim the benefit under Title 35. United States 
Code. §120 of any United States application (s) listed 
below and, insofar as the subject matter of each of the 
claims of this application is not disclosed in the prior 
United States application in the manner provided by 
the first paragraph of Title 35, United States Code, 
§122, I acknowledge the duty to disclose material 
information as defined in Title 37, Code of Federal 
Regulations, § 1.56(a) which occured between the filing 
date of the prior application and the national or PCT 
international filing date of this application. 



PCT/DE00/02434 
(Application Serial No.) 
(Anmeldeseriennummer) 



19.07.2000 

(Filing Date D, M, Y) 
(Anmeldedatum T, M, J) 



(Status) 

(patentiert, anhangig, 
aufgegeben) 



(Status) 

(patentiert, anhangig, 
aufgeben) 



pending 

(Status) 

(patented, pending, 
abandoned) 



(Status) 

(patented, pending, 
abandoned) 



Ich erklare hiermit, dass alle von mir in der vorliegen- 
den Erklarung gemachten Angaben nach meinem 
besten Wissen und Gewissen der vollen Wahrheit 
entsprechen, und dass ich diese eidesstattliche Erkla- 
rung in Kenntnis dessen abgebe, dass wissentlich und 
vorsatzlich falsche Angaben gemass Paragraph 1001, 
Absatz 18 der Zivilprozessordnung der Vereinigten 
Staaten von Amerika mit Geldstrafe belegt und/oder 
Gefangnis bestraft werden koennen, und dass derartig 
wissentlich und vorsatzlich falsche Angaben die Gul- 
tigkeit der vorliegenden Patentanmeldung oder eines 
darauf erteilten Patentes gefahrden konnen. 



I hereby declare that all statements made herein of my 
own knowledge are true and that all statements made 
on information and belief are believed to be true, and 
further that these statements were made with the 
knowledge that willful false statements and the like so 
made are punishable by fine or imprisonment, or both, 
under Section 1001 of Title 18 of the United States 
Code and that such willful false statements may 
jeopardize the validity of the application or any patent 
issued thereon. 
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Patent and Trademark Office-U.S. DEPARTMENT OF COMMERCE 



German Language Declaration 



POWER OF ATTORNEY: As a named inventor, I 
hereby appoint the following attomey(s) and/or 
agent(s) to prosecute this application and transact all 
business in the Patent and Trademark Office 
connected therewith, (list name and registration 
number) 

■nun 

29172*, hereby appoint 

progress patent trademark office 



Direct Telephone Calls to: (name and telephone 
number) 

Ext. 



Postanschrift: Send Correspondence to: 

Bell, Boyd & Lloyd LLC 

Three First National Plaza, 70 West Madison Street, Suite 3300 60602-4207 Chicago, Illinois 
Telephone: (001) 312 372 11 21 and Facsimile (001) 312 372 20 98 
or 

Customer No. in progress 



Voller Name des einzigen oder urspriinglichen Erfinders: 

Dr. ANDREAS JURISCH i - 


Full name of sole or first inventor: 

Dr. ANDREAS JURISCH 


Unterschrift des Erfinders Datum 


Inventor's signature Date 


Wohnsitz / 

SCHWANTE, Deutschland 


Residence f 

SCHWANTE . GE RMAN Y X><s^ 


Staatsangehorigkeit 

DEUTSCH 


Citizenship 

GERMAN 


Postanschrift 

EICHENWEG 11 


Post Office Addess 

EICHENWEG 11 


16727 SCHWANTE 
DEUTSCHLAND 


16727 SCHWANTE 
GERMANY 


Voller Name des zweiten Miterfinders (falls zutreffend): 

Dr. DIETER KRAMER 2 - «o© 


Full name of second joint inventor, if any: 

Dr. DIETER KRAMER 


Unterschriftdes-fefinders / / Datum 


Seconer1nvento£s signatuoi/' Date 


DRESDEN, DEUTSCHLAND 


DRESDEN.GERMANY i3€L/0 


Staatsangehorigkeit 

DEUTSCH 


Citizenship 

GERMAN 


Postanschrift 

RITTERSHAUSSTR. 18 


Post Office Address 

RITTERSHAUSSTR. 18 


01279 DRESDEN 
DEUTSCHLAND 


01279 DRESDEN 
GERMANY 



(Bitte entsprechende Informationen und Unterschriften im (Supply similar information and signature for third and 
Falle von dritten und weiteren Miterfindern angeben). subsequent joint inventors). 



VERTRETUNGSVOLLMACHT: Als benannter Erfinder 
beauftrage ich hiermit den nachstehend benannten 
Patentanwalt (oder die nachstehend benannten 
Patentanwalte) und/oder Patent-Agenten mit der 
Verfolgung der vorliegenden Patentanmeldung sowie 
mit der Abwicklung aller damit verbundenen Geschafte 
vor dem Patent- und Warenzeichenamt: (Name und 
Registrationsnummer anfuhren) 



Telefongesprache bitte richten an: 
(Name und Telefonnummer) 
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